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The cytoplasm and biomembranes in biological cells contain large numbers of proteins that cycli¬ 
cally change their shapes. They are molecular machines that can function as molecular motors or 
carry out many other tasks in the cell. We analyze the effects that hydrodynamic flows induced 
by active proteins have on other passive molecules in solution or membranes. We show that the 
diffusion constants of passive particles are enhanced substantially. Furthermore, when gradients of 
active proteins are present, a chemotaxis-like drift of passive particles takes place. In lipid bilayers, 
the effects are strongly nonlocal, so that active inclusions in the membrane contribute to diffusion 
enhancement and the drift. The results indicate that the transport properties of passive particles in 
systems containing active proteins machines operating under nonequilibrium conditions differ from 
their counterparts in systems at thermal equilibrium. 


I. INTRODUCTION 

Protein machines play a fundamental role in biologi¬ 
cal cells Operating as motors, they are responsible for 
intracellular transport and force generation. As manip¬ 
ulators, they perform various operations involving other 
biomolecules, including RNA and DNA. As pumps, they 
transfer ions across biomembranes. A common feature 
of protein machines is that they undergo cyclic confor¬ 
mational changes that are induced by ligand binding 
and product release. Thus, protein machines are en¬ 
zymes where substrate binding, catalytic conversion to 
products and product release are accompanied by inter¬ 
nal mechanochemical motions. Conformational changes 
within turnover cycles are also characteristic of many 
other enzymes, which need not function as molecular ma¬ 
chines. The results we present in this paper are also ap¬ 
plicable to such enzymes. 

When a macromolecule cyclically changes its shape, it 
induces hydrodynamic flows in the surrounding fluid or 
biomembrane in which it resides. Such pulsating flows 
can act on other passive particles in solution or lipid bi¬ 
layers. The aim of the present study is to analyze the col¬ 
lective hydrodynamic effects that active macromolecules 
have on passive particles in the medium. We shall show 
that these effects lead to substantial modifications of the 
diffusion constants of passive particles. Furthermore, di¬ 
rected drift of passive particles can be induced when there 
are spatial gradients of active macromolecules, a phe¬ 
nomenon that is reminiscent of chemotaxis. 

The investigation of hydrodynamic effects in active flu¬ 
ids is an important field of current research.® 4 While the 
hydrodynamics of bacterial motion has been studied of¬ 
ten, the active elements may be of inorganic origin and 
operate through various flow-generation mechanisms.® ' 13 ' 
There has been a considerable amount of work on swim¬ 
mers that can propel themselves by cyclically changing 


their shapes.^ Interactions between such swimmers and 
their collective flows have been analyzed DSH2D Also, inter¬ 
actions between active hydrodynamic dipoles ha ve bee n 
investigated t heoret icall y 2 2 * 23 1 and experimentallyP^EI! 

Hydrodynamic effects on individual protein machines 
have been studied; for example, investigations of simple 
models of such machines have shown how they propel 
themselves and behave under a loacP^ also, the effects 
of hydrodynamic interactions on the internal dynamics 
have been analyzed 2 -. Active protein inclusions in lipid 
bilayers can act as hydrodynamic dipoled 2 ® and, under 
certain conditions, such inclusions can behave as active 
membrane swimmers 29 . 

The focus of the present study differs in several re¬ 
spects from the work recounted above. We are not be 
interested in the effects of hydrodynamics on the opera¬ 
tion of a single machine, nor do we consider interactions 
between protein machines and the collective flows exhib¬ 
ited by such hydrodynamically-coupled active molecules. 
Instead, we concentrate on the effects that protein ma¬ 
chines can have on passive particles in the medium. 
While some proteins can indeed behave as motors, we 
only require that such machines act as cyclic hydrody¬ 
namic dipoles. Consequently, our analysis concerns the 
statistical effects that populations of incoherently oscil¬ 
lating dipoles can have on passive particles in the system. 

In bulk solution, the laws of three-dimensional (3D) 
hydrodynamics need to be applied; however, as already 
pointed out by Saffmann and DelbriidP^, biological 
membranes should behave as two-dimensional (2D) flu¬ 
ids when lipid flows in a membrane that occur on scales 
shorter than a micrometer are considered. 30 32 Recently, 
2D lipid flows were directly observed in mesoscopic sim¬ 
ulations of lipid bilayers^S additionally, 2D diffusion in 
biomenrbranes was experimentally investigated 3 ^. It is 
well known that 2D hydrodynamics is characterized by 
the presence of ultra-long-ranged logarithmic interactions 






2 


that make it qualitatively different from the 3D easeful 
We study the hydrodynamic effects of active machines in 
both 3D and 2D systems. 

The outline of the paper is as follows. First, we present 
a simple model of a protein machine as a dimer under¬ 
going random cyclic changes. This model is only used 
as an illustration and our general results do not rely on 
this specific structure of the model. We then show that 
a cyclically fluctuating hydrodynamic dipole will induce 
diffusive motion and directed drift of a passive particle 
located at some distance from it. When a population 
of cyclic hydrodynamic dipoles is randomly distributed 
in medium, they will enhance the diffusion of all pas¬ 
sive particles in the medium. Moreover, if such dipoles 
are non-uniformly distributed and concentration gradi¬ 
ents in these species are present, directed flows of all 
passive particles will be induced. Numerical estimates of 
the magnitudes of the effects are given, and a discussion 
of the results is presented. 


II. PROTEIN MACHINES AS FORCE DIPOLES 

Molecular machines are biomolecules, most often pro¬ 
teins, that undergo structural changes in shape during 
their operation cycles. These shape changes, induced 
by ligand binding and product release, take place un¬ 
der nonequilibrium conditions; therefore, they differ from 
thermally-induced shape fluctuations for which micro¬ 
scopic reversibility holds and the fluctuation-dissipation 
theorem applies. These molecular machines operate in a 
viscous environment and their dynamics takes place un¬ 
der low Reynolds number conditions so that inertia does 
not play a significant role. As a result, if a force is ap¬ 
plied to a particle in the fluid, the same force acts on 
the fluid. More specifically, such protein machines act as 
stochastic oscillating force dipoles that can influence the 
motions of other particles in the system. 

In order to illustrate how oscillating force dipoles arise 
in protein machines, consider a coarse-grain model of a 
protein comprising two domains, represented by beads, 
and connected by an elastic spring. If the spring con¬ 
tracts or stretches, two forces equal in magnitude and op¬ 
posite in the direction will act on the beads, and also on 
the surrounding fluid, creating a force dipole of strength 
m = FI , where F is the magnitude of the force and l 
is the distance between the beads. As a consequence of 
ligand binding and product release the characteristics of 
the spring force will change. 

A simple model for such ligand-induced force changes is 
as follows. The overdamped motion of the beads, labeled 
by 1 and 2, is given by ±\ = pk(— X\ — Z(s)), where 
fi is the mobility of the beads, k is the stiffness constant 
of the elastic spring, and l(s) is the natural spring length 
which depends on the variable s = {0,1} that specifies 
the binding state of the model machine: l(s = 0) = a and 
l(s = 1) = b . If the binding state of the machine is fixed, 
it reaches a steady configuration with the spring length 


a or b and remains indefinitely in that state. Suppose 
however that the internal state can change as a result 
of reactions, so that stochastic transitions from s = 0 
to s = 1 and back take place. Generally, such transi¬ 
tions will depend on the length of the spring: close to 
the state with l = a the transition from s = 0 to s = 1 
becomes possible, whereas near the state with l = b the 
reverse transition can occur. Under these conditions, the 
machine will go through cycles of expansion and contrac¬ 
tion, similar to the cyclic conformational changes in a real 
protein machine. It is an oscillating force dipole which 
can act on the surrounding viscous fluid to create hydro- 
dynamic flows that can induce motions in other passive 
particles 


III. HYDRODYNAMIC EFFECTS 


When a force F is applied to the fluid at a point r, it 
produces a fluid flow field at R that advects a particle at 
this location with velocity, 

R a = G aP (R-T)Fp( r ), (1) 


where G a @ is the mobility tensor which, for sufficiently 
large distances, can be evaluated in the Oseen approx¬ 
imation. (The Einstein summation convention will be 
used throughout this paper.) For an oscillating dimer of 
length x with orientation given by the unit vector e and 
interaction force F, we have R a = [G a/ ?(R — r — xe) — 
G ai g(R — r)\epF. If the dimer length x is small compared 
to the distance |R — r'| , we can write 


dG af) 

— e,e 7 m, 


( 2 ) 


where m(t) = x(t)F(t) is the force dipole. 

Although we derived eq. ([2]) for a specific model, it is 
general. If an object immersed in the fluid changes its 
shape, it generates a hydrodynamic flow that, at large 
separations from this object, can be described as be¬ 
ing produced by a force dipole. Below, we shall treat 
a protein machine as a stochastic force dipole m(t) with 
zero mean, (m(t)) n = 0, and correlation function S(t) = 
(m(t)m(0))„, where (•••)„ denotes an average over the 
stochastic fluctuations of the force dipole. According to 
eq. a force dipole will induce motion in a passive 
particle in the fluid. 

Consider a collection of such force dipoles in the sys¬ 
tem, located at positions {RJ with orientations {e.;}. If, 
within the time interval being considered, the displace¬ 
ments in the position R(f) of a passive particle are small, 
we can write R a (t) = Ro, a + Pa(t) and retain only terms 
that are first order in the displacements p(t ) , so that 


Pa — 




dG a0 ( Ro-Ri) , d 2 G a p(R, 0 - Ri) 


dRj 

L ’ 

xei,pe i:1 mi(t). 


dRi^dRig 


PS 


( 3 ) 
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It is convenient to write this equation in field-point no¬ 
tation, 


Pa 


Y- f , raG„ s (R 0 -r) a 2 G„j(R«-r) i 

2 r J dr l - Bfr -+ — Syg7 s —«J 

X &i ,/?e^ 7 77q (t)<5(Rj I*), (4) 

i 


where the last term represents the microscopic density of 
the product e^e^TO^t) at r. Note that the first and the 
second derivatives of the Green function correspond to 
dipole and quadrupole terms. This equation shows that 
the instantaneous position of a passive particle evolves 
with time according to a stochastic differential equation 
with multiplicative noise. We assume that the fluctu¬ 
ations of the different force dipoles are uncorrelated so 
that (mi(t)rrij( 0)) n = S(t)Sij. 

The diffusion constant D and mean drift velocity Vof 
the passive particle under the action of a field of fluctu¬ 
ating hydrodynamic dipoles can be determined from 


D= ^l V a = (V a ), (5) 

where d is the dimensionality of the system, SV a = V a — 
V a and the angle bracket (• • •) denotes an average over 
the stochastic fluctuations of the force dipoles, as well as 
their orientations and positions. 

The expression in eq. ([I]) for the velocity V a = p a of 
a particle may then be substituted into these equations, 
retaining only the leading terms, to obtain the diffusion 
coefficient and drift velocity. When computing the aver¬ 
age values in eqs. ([5]) we assume that the orientations of 
the force dipoles are not correlated with their positions 
so that epep>e^ey5(Ri — r)) = flpp^yc^), where 
f Ipp^y = (epep>e-yey) and c(r) = (S(R,; — r)) is the 

local concentration of force dipoles at point r in the fluid. 
We find 


D(Rq) = ^^/3/3'7 7 ' (6) 

dG aP { Rp-r) 0G a fl,(Ro-r) 


dr 


dry 


V q (R 0 ) = Sfiflfj/. 


dr 


d 2 G a p{ Ro-r) dG^Ro-r) 


dr-ydrs 


dry 


(7) 

:( r ). 


Here S = / 0 °° dt S(t). We shall now analyse these equa¬ 
tions separately for 3D and 2D systems. 


( 8nr/r)- 1 (5 a /3 + f a fp) where g is the fluid viscosity and r 
is the unit vector specifying the direction of r. Suppose 
that the machines are uniformly distributed in space and 
their average concentration is cq. Then, eq. § yields 


Sc 0n [' dG a p(r) dG a p>(r) 
= —Q d0/j y J dr— - - - 


7' 


( 8 ) 


In writing this expression we have made the change of 
integration variable, R 0 — r —> r, retaining the same sym¬ 
bol to avoid proliferation of notation. The prime on the 
integral indicates that the integral is cut off at short dis¬ 
tances, since overlap between the passive particle and an 
active enzyme is restricted by volume exclusion. We de¬ 
note this cut-off distance by £ c and choose it to be larger 
than the mean diameter of the passive and active parti¬ 
cles to account for the fact that the Oseen approximation 
breaks down at distances larger than the particle sizes. 
Introducing dimensionless coordinates z = r/£ c , the dif¬ 
fusion coefficient can be estimated as 


where the positive factor £ is given by 

>_ ^/3/3'77' f , dg a p(z) dg af) > (z) 

1927r 2 J dz 1 dzy 


(9) 


( 10 ) 


and the dimensionless function g a p( z) = z^ 1 {8 a p-\-z a zp). 

Turning to eq. ([7]), we notice that the drift velocity 
vanishes for a uniform distribution of active particles. 
Suppose instead that a constant concentration gradient 
in the direction n, Vc = n(n ■ Vc) = n(Vc), is present 
and c(r) = Cq + (r ■ n)Vc. Now, the integration in eq. ([7]) 
yields 




dr 


d 2 G a0 ( r) dGsp'ir) 

dr^drs dry 


(r-n)Vc. (11) 


Again, introducing dimensionless coordinates z = r/£ c . 
we find 


V 



£Vc, 


where the positive dimensionless factor £ is 


( 12 ) 


£ 


^I/3/3 / 77 / 

647t 2 


dz n a 


d 2 g ai g(z) dgs/3' (z) 

dz-ydzs dzy 


(z • n). 


(13) 


B. Two-dimensional systems 


A. Three-dimensional systems 

For applications to protein machines in bulk solution, 
for example in the cytoplasm of biological cells, the 3D 
Green function in the Oseen approximation is G a p{ r) = 


As noted earlier by Saffman and DclbriiclP^, biolog¬ 
ical membranes should behave as two-dimensional lipid 
fluids on sub-micrometer length scales. Therefore, the 
effects of an ensemble of active protein machines on the 
motion of a passive particle in a lipid bilayer is an exam¬ 
ple where a 2D description is appropriate. In the Oseen 
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approximation, the 2D Green function of lipid bilayers 


G af g(r) = (47r7y m ) a (-(l + In (hir))5 a p + r a rp) (14) 

Here, ?y m is the two-dimensional viscosity of the lipid bi¬ 
layer, which is related to its 3D viscosity rjp by r] m = ly^h, 
where h is the thickness of the bilayer. In contrast to the 
3D case, hydrodynamic interactions in 2D are ultra-long- 
ranged, due to the logarithmic dependence on the dis¬ 
tance r. For biomembranes, one can use the estimate 31 
k _1 = r)ph/{2if) where rj is the viscosity of the surround¬ 
ing aqueous medium. Typically, k” 1 is of the order of a 
micrometer. At separations larger than this distance, 3D 
effects become important for biomembranes. Our anal¬ 
ysis applies only for relatively small membranes of mi¬ 
crometer size; for larger membranes the crossover to 3D 
hydrodynamics at long length scales needs to be taken 
into account. 

The general expressions in eqs. ([6]) and Q hold in the 
2D case as well; however, because of the logarithmic dis¬ 
tance dependence in the Green function, hydrodynamic 
effects are nonlocal. Therefore, it is not possible to ob¬ 
tain estimates similar to those in eqs. ([£]) and (fl2| for 


such systems. Still, some general estimates can be made. 

Suppose that the concentration distribution of active 
inclusions has the form c(r) = Cox(r) where x( r ) is of or¬ 
der unity and varies on some characteristic length scale 
£0 ~ k - 1 . Substituting this into eq. <[tJ) and introduc¬ 
ing the rescaled dimensionless coordinates z = r/£ 0 , we 
obtain 


D( R 0 ) « ^C 2 (Ro), 

where the dimensionless function £ 2 (Ro) is given by 

C 2 ( 0 ) 327r 2 


(15) 


(16) 


x j"' dz dq a p{ Rq/Iq - z) dq a p>( Rq/Iq - z) 


dz~ 


dzy 


and q af 3 (z) = - In {z)6 a p + z a zp. 

In a similar manner, the drift velocity in 2D systems 
can be estimated: 


F a (R 0 ) » /^_C 2 ,a(Ro), 


where 


£,2,a (Rq) = 


Q 


167T 2 


(17) 


(18) 


x I"' dz d 2 q a0 (R o /£ o - z) dqgp'j Rp/4 - z) ^ 


dzrydzs 


dZ-y 


An important difference with the 3D case is that now 
the diffusion constant and the drift velocity are not deter¬ 
mined only by the local concentration of the protein ma¬ 
chine and its gradient. Instead, they depend on the con¬ 
centration distribution over the entire 2D system through 


the dimensionless functions £2 and £ 2 , a- The integrals 
eqs. (161 and (18) diverge logarithmically at large dis¬ 


tances and thus the finite size of the membrane is impor¬ 
tant. 


IV. NUMERICAL ESTIMATES 

The magnitude of a force dipole m of a protein ma¬ 
chine can be roughly estimated as m ~ Fdp , where F is 
the force generated by the machine and dp is the linear 
size of the protein. Molecular motors, such as myosin 
or kinesin, typically generate forces about 1 pN and this 
can be chosen as the characteristic value for F. Taking 
the size of a protein to be about 10 nm, the force dipole 
can be estimated to be m = 10~ 20 N m. The correlation 
time for force-dipole fluctuations can be taken to be the 
duration t c of the turnover cycle in a chemical machine. 
We choose a time of about t c ~ 10 ms. The parame¬ 
ter S can then be evaluated to give S ~ m 2 t c = 10 -42 
N 2 m 2 s. Note that this estimate corresponds to substrate 
(ATP) saturation conditions: the machine binds a new 
substrate molecule and enters into a new cycle immedi¬ 
ately once the previous cycle finishes. If this condition 
is not satisfied, the protein machine must wait for a new 
ATP molecule to arrive. During this waiting period, the 
machine does not act as a force dipole and this will de¬ 
crease the value of S. 

Concentrations of active proteins inside a biological cell 
can vary over a large range. The highest concentrations 
of the order of 10 -4 M are characteristic for the enzymes 
involved in glycolysis. As a rough estimate, a value of 
1(V 6 M can be chosen, corresponding to Cq = 10 21 m -3 , 
so that the mean distance between the proteins is about 
100 micrometers. Given the protein size, we choose a 
cutoff length of £ c = 20 nm. The viscosity of water is 
about 10 -3 Pa s. The dimensionless numerical factors in 


eqs. (101 and (13) can be computed and take the approx¬ 
imate values £ = 1CT 2 and £ = 10 -3 , 

With these values, the contribution § to the dif¬ 
fusion coefficient due to hydrodynamic effects arising 
from protein machines in bulk 3D solutions is about 
D = 10~ 7 cm 2 /s. This result should be compared with 
typical diffusion constants in water that can vary from 
about 10 _5 cm 2 /s for small molecules to 10~ 7 cm 2 /s for 
small proteins of the size considered here in water. To 


estimate the magnitude of the drift velocity from eq. (12), 


we can take Vc = A c/L where Ac is the concentration 
difference across the cell and L = 10 /im is the typi¬ 
cal eukaryotic cell size. If we take Ac ss O.lco, where 
again Co = 10 21 m -3 , we obtain a drift velocity of about 
V a = 1 qm/s. 

Proceeding to lipid bilayers, we observe that their 3D 
viscosity t)l ~ 1 Pa s is typically a factor of 10 3 higher 
than that of water. The 2D viscosity of such bilayers is 
Vm = r/ph where h ~ 1 nm has been taken as the thick¬ 
ness of the bilayer. The magnitude S ~ 10 -42 N 2 m 2 s 
may again be used for the force dipoles. Taking the 
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mean distance between proteins to be approximately 100 
nm, the 2D concentration is about Co = 10 14 m~ 2 . The 
dimensionless factors £2 and £ 2 ,a depend on the con¬ 
centration distribution c(r) of active protein inclusions 
over the entire membrane, according to eqs. (161 and 
(18). As rough estimates of these factors we choose 
C 2 = 6,a = IQ" 2 - 


Given these numerical values, the hydrodynamic ef¬ 
fects of active protein inclusions are predicted to increase 
the diffusion of passive particles within the membrane by 
about D = 10~ 8 cm 2 /s . For comparison, thermal diffu¬ 
sion constants for proteins in lipid bilayers are of the 
order of 10~ 10 cm 2 /s and diffusion constants for lipids 
are about 10~ 8 cm 2 /s. The characteristic drift velocity of 
passive particles in lipid bilayers is estimated to be about 
V a = 1 /rm/s , assuming that the characteristic length 
for concentration variation in a membrane is about one 
micrometer. 

These numerical estimates should be used only as 
rough guide to the possible magnitudes of the effects 
because many of the parameters may vary significantly 
from one system to another, or are only known poorly. 
For example, forces have only been measured for some 
molecular motors and, for protein machines that are not 
motors or for enzymes, they may be smaller than one 
pN. However, the concentrations of some proteins that 
behave as force dipoles may be significantly higher than 
the value we have assumed. For example, the enzyme 
phosphoglycerate kinase involved in glycolysis is present 
in the living cell in the concentrations up to 10~ 4 M, 
two orders of magnitude larger than our assumed value. 
Proteins typically contribute about 40% of the mass in 
biomenrbranes and their concentrations may well be sig¬ 
nificantly higher than our 2D values. There is also some 
ambiguity in the choice of the cutoff length in 3D sys¬ 
tems but its lower limit is controlled by the distance cor¬ 
responding to the volume exclusion of passive and active 
particles. Consequently, the uncertainty in the estimates 
of hydrodynamic effects is high and deviations of up to 
two orders of magnitude from the numerical estimates 
given above may well be possible. Moreover, these es¬ 
timates hold only under conditions where the substrate 
(ATP) is in excess. If the rates of machine cycling are 
limited by substrate supply, the effects are weaker and 
obviously vanish in absence of the substrate. 


V. DISCUSSION AND CONCLUSIONS 

In a biological cell there are large populations of ac¬ 
tive proteins, both molecular motors and enzymes that 
change their conformations within catalytic cycles. In 
this paper we showed that when active proteins are 
present, either in solution or in lipid bilayers, they can 
substantially modify the diffusion constants of passive 
particles in the system. These modifications affect all 
passive particles, and all active proteins, even of different 
kinds, contribute to the effect provided they are supplied 


with substrate and remain active. The magnitude of the 
effect can be comparable to the value of thermal diffu¬ 
sion constants under physiological conditions. Further¬ 
more, if protein machines are nonuniformly distributed 
in a cell or in a biomembrane, directed drift of passive 
particles, analogous to chemotaxis, can occur; however, 
the mechanism is completely different: all active proteins 
contribute towards it and all passive particles experience 
the drift. Drift velocities of the order of micrometers per 
second can be realized. The enhancement of diffusion 
and cheinotaxis-like drift should take place for the pro¬ 
tein machines (enzymes) themselves as well and we note 
that enhanced diffusion of catalytically-active enzymes 
has bee n obs erved and is a topic of considerable current 
interest E5122I 

In three dimensions, hydrodynamic interactions are al¬ 
ready long-ranged, with power-law dependence. In two 
dimensions they become ultra-long ranged with a loga¬ 
rithmic dependence on distance. Thus, the effects pre¬ 
dicted to exist in 3D and 2D systems differ substantially. 
In solution, the change in the diffusion constant is de¬ 
termined by the local protein machine concentration and 
the drift velocity is controlled by its spatial gradient. In 
contrast, in 2D systems such as lipid bilayers, the effects 
are essentially nonlocal: the change in the diffusion con¬ 
stant and the drift are determined by the distribution of 
active inclusions over the entire membrane. 

Our description of how active molecular machines, 
through hydrodynamic interactions, influence the dy¬ 
namics of passive particles was based on the equations 
of continuum hydrodynamics. One might question the 
use of such a continuum description for molecular sys¬ 
tems. It is well established 37 that hydrodynamic effects 
are observable on very small scales of tens of solvent par¬ 
ticle diameters. They persist in spite of strong fluctua¬ 
tions and their presence is signalled in the long-time tails 
of velocity correlation functions or even in the transport 
properties of polymers. Our use of continuum equations 
is restricted to rather long scales so that the main con¬ 
clusions of our study should be robust. 

Our analysis has focused on general aspects of the phe¬ 
nomena and is not intended to address the full complex¬ 
ity seen in biological systems. The cell is crowded with 
obstacles of various kinds and these crowding elements 
could modify the hydrodynamic behavior. Nevertheless, 
the results suggest a physical picture of the dynamics 
in the cell that differs from one based solely on thermal 
motions. Whenever active proteins are present, they not 
only perform their varied biological functions, but also 
generate hydrodynamic fields that span the entire cell. 
Such non-thermal fluctuating fields affect the diffusion of 
particles and lead to their directed drift. There is an¬ 
other important effect, which we have not analyzed, but 
is worth noting: because such fluctuating fields repre¬ 
sent non-thermal noise, they can be rectified and work 
or energy can be extracted from them. Through hydro- 
dynamic interactions, active inclusions supply power in a 
distributed way within the cell or a biological membrane. 
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